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Abstract 

This paper presents an event-driven way finding algorithm for pedestrians in an evacuation scenario, which 
operates on a graph-based structure. The motivation of each pedestrian is to leave the facility. The events 
used to redirect pedestrians include the identification of a jam situation and/or identification of a better route 
than the current. This study considers two types of pedestrians: familiar and unfamiliar with the facility. 
Four strategies are modelled to cover those groups. The modelled strategies are the shortest path (local and 
global); They are combined with a quickest path approach, which is based on an observation principle. In the 
quickest path approach, pedestrians take their decisions based on the observed environment and are routed 
dynamically in the network using an appropriate cost benefit analysis function. The dynamic modelling of 
route choice with different strategies and types of pedestrians considers the manifold of influences which 
appears in the real system and raises questions about the criticality of an evacuation process. To address 
this question criteria are elaborated. The criteria we focus on in this contribution are the evacuation time, 
the individual times spent in jam, the jam size evolution and the overall jam size itself. The influences 
of the different strategies on those evaluation criteria are investigated. The sensibility of the system to 
disturbances (e.g. broken escape route) is also analysed. 
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1. Introduction 

In recent years pedestrian dynamics has gained more importance and a lot of attention due to con- 
tinuously growing urban population and cities combined with an increase of mass events. This sets new 
challenges to architects, urban planners and organizers of mass events. One of the main goals is an effective 
use of the designed facility, for instance by minimizing jams thereby optimizing the traffic flow. In this con- 
text pedestrian simulations are already used, e.g. for escape route design |Tl[2j[3l[4]. However the approaches 
used in these simulations are only the first step to model the manifold influences on human beings during 
an evacuation process. Major issues in this area include orientation and way finding: Given a set of possible 
routes, which criteria influence pedestrians choice for a particular route? This is essential for reproducing 
route choice in computer models and is difficult due to the many underlying subjective influences on this 
choice. The manner by which pedestrians choose their way has a direct influence not only on the overall 
evacuation time but also on the average time pedestrians spend in a jam. In this paper we restrict ourselves 
to the case where pedestrians all have the same motivation, to leave the facility. 

The approaches of modelling pedestrians motion fall into two main groups: microscopic and macroscopic 
models. Microscopic models are further categorized in spatially discrete (e.g. Cellular Automata [5l[6l[7]) 
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and spatially continuous models (e.g. Social Force Model [8 , Generalized Centrifugal Force Model [9]). For a 
detailed overview we refer to [10 . CA models use floor fields (static and/or dynamic) to direct pedestrians to 
a destination point. Route choice in continuous models can be achieved by means of a network which consists 
of a set of destination points. This type of way finding is known as graph-based routing p!T| \T2\ {T3\ [M) [T5]. 
The destination points can be pre-determined (exits for instance) or adjustable (crossings, turning point at 
the end of a corridor for instance). The minimal network is usually a a visibility graph (see [16 for more 
details) which ensures that any location on the facility is within the visibility range of at least one node. 
The initial network can be refined by adding more adjustable points converging to CA [17]. This graph can 
be extended to a navigation graph by adding more adjustable points. The generation of such graph is a 
complex process, some efficient methods are presented in [18l[T9]. Once built the shortest path is usually 
determined using well established algorithms such as Dijkstra [20^ or Floyd warshall [21]. Such networks are 
widely spread in motion planning by robots as well. 

The intrinsic behaviour of humans in the case of an evacuation is generally to follow the seemingly 
(self estimated) quickest path. This is indeed a subjective notion as it depends on some prerequisites, e.g. 
whether or not the pedestrian is familiar with the facility. The modelling of the quickest path is achieved 
by systematically avoiding congestions. In CA this is achieved by means of dynamic floor fields [22] where 
pedestrians moving increase the probability of using that path thereby making it more attractive for other 
pedestrians. This implicitly leads to congestion avoidance. The density in front of the moving pedestrians 
within their sight range is also considered [23] as well as the payload at exits [24 ; other approaches include 
navigation fields [25^, ^6 . Continuous models usually optimize the travel time in the constructed network. 
In [27^, pedestrians minimize their travel time by solving the Hamilton- Jacobi-Bellman equation yielding 
to the optimal pedestrians path at each time step. A combination of a graph-based routing with CA is 
presented in [18 where the fastest path for pedestrians is computed using a heuristic A* algorithm [28 . 

In this contribution an event driven way finding in a graph-based structure is introduced. The approach 
is based on an observation principle, pedestrians observe their environment and take their final decision 
based on the obtained data. With this observation the quickest path is achieved. The modelled strategies 
are the local shortest path (LSP), the global shortest path (GSP), the local shortest combined with the 
quickest path (LSQ) and finally the global shortest combined with the quickest path (GSQ). 

An important point less discussed is given by the criticality of an evacuation process, first and foremost 
the meaning of criticality for an evacuation process. Evaluation criteria like the building itself, the population 
size and the initial distribution of the evacuees, the evacuation time are discussed in [29]. More individual 
criteria like the individual travel time and waiting time are investigated in [30 . The most used criteria are 
the overall evacuation time and a visualisation of the evacuation process at specific times. In this paper we 
elaborate other criteria to address the criticality of an evacuation simulation. The analysed criteria are the 
individual time spent in jam, the jam size evolution over time and the total jam size defined as the area 
under the jam size evolution. We give a special credit to the time pedestrians spend in jam as well as the 
jam size. This is of particular importance. In the case of the Hermes project [31], there are congestions that 
vary (in place and size) depending on the type of people attending the events. We try to reproduce this 
observed phenomenon for the forecast of the evacuation dynamic. In order to achieve this, route choice has 
to be individually modelled. Those events involve pedestrians familiar and unfamiliar with the facility. The 
four previously mentioned strategies are proposed to reproduce their route choice. Their influences on the 
previous mentioned criteria are investigated. 

The modelling approaches are presented in the second section. The third section describes the evacuation 
assessment criteria. The analysis of the results including simulation, distribution of the evacuation time and 
distribution of the individual times spent in jam using different initial conditions are discussed in the fourth 
section. 

2. Modelling 

The framework used for describing pedestrian traffic can be divided in a three-tiers structure. One 
distinguishes between the strategic, the tactical and the operational level [27]. In our model the operational 
level of the pedestrian walking is described by the Generalized Centrifugal Force Model [9] which operates 
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in continuous space. In this model pedestrians are described by ellipses. The semi-axes of the ellipses are 
velocity dependent, faster ellipses need more space in the moving direction. The fundamental diagram is 
reproduced by the model at corridors making it adequate for the analysis presented here. In this paper we 
focus on the strategic level only, i.e. the pedestrians are solely given the next destination point, which is the 
next intermediate destination for the self estimated optimal route. This is also the direction of the driven 
force. We simulate pedestrians that are familiar with the facility and pedestrians unfamiliar with the facility. 
For those two groups we identify criteria to reproduce their route choice. Those criteria are the local/global 
shortest path and the quickest path. The dynamic change of the strategies is modelled. This emulates the 
internal state of the pedestrians and the strategies are subject to change during the evacuation process. In 
addition one of the challenges consists of finding a good balance between the number of parameters and 
the numbers of criteria. The model should be as simple as possible with the parameter space kept as small 
as possible while considering as many criteria as required. This is important for model understanding and 
stability. 

In the framework used here, pedestrians move from one decision area to the next one. The decision 
areas are connected with nodes, which will be interchangeable with destination points. A decision area is 
a place where the pedestrian decides which way to go or change their current destination. In this work a 
decision area is an abstraction for rooms and corridors, in addition we restrict ourselves to the case where 
the destination points are exits and corridors end. Fig. [l] illustrates this principle and is a mapping of the 
facility presented in Fig. [2] The pedestrians will be moving from the decision area 1 to the decision area 2. 
The two areas are connected with one node rii. The network is automatically generated from the facility 
based on the inter-visibility of the exits. The Euclidean distance between the nodes is used as weights. 

2.1. Shortest Path 

The most straight-forward routing approach is the local shortest path. Once a node in the network 
is reached, the local nearest node is chosen as next destination. The global shortest path is determined 
by running the Floyd Warshall algorithm with path reconstruction [21 on the built graph network. The 
runtime of (9(n^) is not an issue since we only have a small amount of nodes. From every node on the graph, 
the global shortest path to reach the outside can be determined. Pedestrians familiar with the facility have 
the global map i.e. the constructed network and may approximate the global shortest path to their final 
destination independently of their current location. They have a better analysis possibility of the current 
situation. Other pedestrians without any global information choose the local shortest path. 

2.2. Quickest Path 

In contrast to the shortest path, the quickest path is dynamic and changes with time throughout the 
simulation. The business logic of the routing algorithm is shown in Fig. |3] The main events used in 
this routing algorithm to redirect pedestrians are the entering a new room and the identification of a jam 
situation. The pedestrians are first routed using the shortest path, global or local depending on theirs 
affiliations [S^- The key elements of a quickest path routing approach are the estimation of the travel time, 
the estimation of the gain and an assessment of this gain. Three functions are developped in this paper to 
model those key elements. We first define four values which will be used throughout this section. 

Definition 2.1. Re-routing time t^ 

The re-routing time tr for a pedestrian i with position Xi is the time where one of the following conditions 
holds: 

\W{t)\\ ^ Vmin^'^t e [tr — tmin^U] V ||:^(tr) " ^ <^min whcrc Vmin ^-^ the threshold jam velocity, train the 
patience time and dmin the minimal distance to the node to consider it as reached. 

Definition 2.2. Reference pedestrian Re f{i^ni) 

Let i and j be two pedestrians with positions Xi and Xj . j is a reference pedestrian to i with respect to the 
node Hi of the graph and is defined as: 

Ref{i, fii) = {j : j e Q{ni)^Xj visible from Xi A\\xj —Xi\\ = min}, where Q{ni) is the jamming queue at the 
node Hi . If more than one pedestrian satisfy the condition, one is randomly chosen. 
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Figure 1: Network mapping of the facility presented in Fig. [2] with 2 decisions areas and 4 nodes. 
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Figure 2: Process of selecting a reference pedestrian prior to a route change. Pedestrians are denoted with their positions, oci 
will select oc^, ocq and ocg- X2 has no clearance of the current situation and will not select any. aJs selects xq and xs- x\ will 
only select xj. 



Definition 2.3. Jamming queue Q{ni) 

Let Hi be a node in the navigation graph, the jamming queue Q{ni) at the node fti is defined as 

Q{ni) = {\/i : Hi destination of the pedestrian i A \\vi\\ < Vmin} where Vmin ^-^ the threshold jam velocity. 

Definition 2.4. Visibility range V{xi) 

V{xi) represents the set of all nodes within the visibility range of the pedestrian i considering the actual 
location Xi in the facility and other pedestrians. It is determined using the Algorithm [7| 



At any time tr (see Definition in 2.1) during the simulation, a new orientation process is started for the 
pedestrian i. A reference pedestrian j (see Definition 2.2) is selected from the queue (see Definition |2.3[ ) and 
observed during an observation time tobs- At the end of the observation, the expected travel time t{xi,nj) 
via all nodes in the visibility range V{xi) (see Definition 2.4) is approximated using Eq. [l] 



t{xi,nj) = < 



11^.11 



pr^j^, if a reference Xj with average velocity Vja was found; 



if the node fij is free (there is no reference in that case). 



(1) 



^obs 



tr~\~t obs 



Vj{t)dt 



(2) 
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Figure 3: Business logic of the quickest path. 
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Vja is the average velocity of the reference pedestrian over the observation time defined by Eq. [2j tobs is 
randomly chosen between 1 and 3 seconds, the minimal distance dmin is set to 0.20 m and the minimal jam 
velocity Vmin is 0.2 ms~^. The estimated travel time is converted to a gain using Eq. [sj 

The cost benefit analysis (cba) function defined in Eq. [4] determines whether it is worth changing the 
current destination, gi and g2 are the gains calculated with Eq. |3] gi is always the current destination 
and g2 the other evaluated alternatives. The benefit returned should be greater than a threshold (gmin) in 
order for the pedestrian to consider the change. The thresholds taken here are 0.20 for familiar and 0.15 for 
unfamiliar pedestrians. 

cba{gi,g2) = ^^—-^ (4) 

When a pedestrian is caught in a jam for a period tmin (patience time) which varies depending on the 
pedestrian, he/she looks for alternatives in the decision area and in the sight range as described in Fig. [2] 
The initial value for tmin is 10 seconds. It is increased by 1 second with any unsuccessful attempt to escape 
the jam. The value is kept until the room is changed. This amortizes the number of routes changes in the 
same room. 

The process of selecting a reference pedestrian is explained in Fig. [2] The pedestrian xi has entered a 
new room. The pedestrians X2, xs and x^ have identified a jam situation. Pedestrian xi selects the reference 
pedestrians X5, xq and xg for the exits n2, ^3 and respectively. The pedestrian X4 has a restricted visibility 
and will have only X7 as reference. There is no possibility for X2 to change. The references selection is based 
on the Euclidean distance and visibility range. The visibility is implemented by drawing a line from the 
concerned pedestrian to the pedestrians in the queue, there should not be any intersections with other 
pedestrians or walls in the room. It is important to mention here, that the queue size does not play a major 
role, more important is the processing speed of the queue. 



Algorithm 1 Visibility Range of Pedestrian Xi'. V{xi) 
Input: Pedestrian Xi 
Output: V{xi) 

V ^ all pedestrians \ {xi} 

W ^ all nodes connected to the actual decision area of Xi 

obstacles JTiin ^ 2 
for all nodes fij in W do 
obstacles ^ 

u ^ queuing pedestrians at node fij 
for all pedestrians Xj in u do 

for all pedestrians in 7^ \ {xj} do 

if [xiXj] intersects with the ellipse shape of the pedestrian Xk then 

obstacles ^ obstacles + 1 
end if 
end for 

if obstacles < obstacles_min and fij ^ V then 

end if 
end for 
end for 
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Unlike other algorithms, the approach presented here is not specialized to a particular case (asymmetric 
exits choice for instance), i.e. not bounded to the geometry. It is also not dependent on the initial distribution 
of the pedestrians. 



3. Evacuation assessment 

Usually evacuation processes are assessed with a visual proof and evacuation time within a feasible range. 
One less discussed question is the criticality of an evacuation process. The state of an evacuation can be 
critical for a certain group of the population, aged persons for example, but rather harmless for a different 
group. Also the same results can be interpreted differently depending on the surrounding conditions. In 
this section we address three criteria to assess an evacuation scenario. 
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(a) Local shortest path 



• 



4m' 



(b) Local shortest with 
quickest path 




(c) Global shortest path 




(d) Global shortest with 
quickest path 



Fi gure 4: Dynamics of the system after 60 seco nds . Co ngestions areas are red. In | (a) | pedestrians follow the local nearest exit. 
In|(c)|the global nearest exit is chosen. In |(b)| and |(d)| pedestrians are allowed to change their current destination. 



3.1. Evacuation time 

The evacuation time is given by the last person leaving the facility. Another definition can be a clearance 
of the building up to 95% of the occupants. We consider the former definition here. A visual assessment of the 
evacuation dynamics for a simple scenario after 60 seconds is given by Fig. |4] The colour of the pedestrians 
is correlated to their current velocity. Red means slow and represents congestions areas. Green corresponds 
to the maximum desired velocity of the pedestrians. The desired velocities are Gaussian distributed with 
mean 1.34 ms"^ and standard deviation 0.26 ms~^ . The evacuation times distribution is presented in Fig. 



5a The smallest variance is achieved by the GSP. This is quite normal as the GSP remains the same for all 



pedestrians independent of their initial positions. The largest variance is given by the LSP. The dynamics 
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brought by the quickest path leads not only to a reduction of the overall evacuation time and a reduction 
of the width but also to realistic shapes in the simulation. 

3.2. Jamming time 

Up to now little importance has been brought to jam analysis itself i.e. how long pedestrian stay in jam. 
This is strongly coupled with the implemented routing strategy. The keyword jam is unfortunately not well 
defined in pedestrian dynamics. A rather crude definition would be to have an absolute zero velocity over a 
minimal time interval. The minimal time interval is needed to avoid very short velocity reduction at sharp 
turn for instance. We consider pedestrians moving at a speed lower than 0.2 ms~^ [33^ for a period of at 
least 10 seconds as being in a jam situation. The total time in jam is recorded for each pedestrian and a 
distribution of the recorded times is calculated. The jamming time distribution for the simulation scenario 
given by Fig. [4] is presented in Fig. [5bj As expected the width of the distribution is smaller for the quickest 
path. 
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Figure 5: Evacuation time and time in jam distribution for 200 pedestrians in the scenario in Fig. ^ 



3.3. Jam size 

The jam size and its evolution cannot directly be derived from the evacuation time and/or the individual 
times spent in jam. It has to be analysed separately and for each exit individually. Short-lived jams are 
rather uncritical, the same holds for short-sized jams. We are particularly interested in big jams with a 
long lifetime. They can be dangerous and should be checked against the characteristics of the evacuees for 
instance their ages, but also against the environmental conditions, the temperature for instance. The jam 
size is calculated by summing up the effective areas occupied by the pedestrians in jam at each time step. 
This is easily done by summing up the ellipses areas representing the pedestrians in the GCFM. Another 
method could be to build the envelope of the pedestrians and calculate the area of the resulting polygon. 
Fig. [6a| shows a straight forward example of a jam size analysis scenario after the first time step (0.01 
second). The colour corresponds to the states. Red pedestrians are waiting for a possibility to move. Other 
pedestrians are already moving. The initial density is on average 1 pedestrian / . The corresponding 
average jam areas are presented in Table [l] The values for the global and local shortest path are the same. 
There as some difference in their corresponding combination with the quickest path. At a first look the 
values for the different routing strategies do not differ much from each other, but one has to consider the 
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difference in the evacuation time, which is more than one minute in this case. There is a distribution of the 
overah load from the congested exits to the non congested ones. 



The jam size evolution at the different exits using the quickest path is shown in Fig. 6b There is a sharp 
increase at the beginning of the simulation. Peaks of approximately 150 are reached at all exits. The 
values then decrease with a slope correlated with the routing strategy used. The slope is constant with the 
static strategies. The quickest path shows a less constant slope due to the fact that some pedestrians are 
changing their destinations. 




(a) Initial homogeneous distribution (b) Jam size evolution at the four exits for the local 

shortest path combined with the quickest path. 



Figure 6: Initial distribution and jam size evolution of 2500 pedestrians in a 50 m x 50 m room with four exits on each wall 
side. The exits width are 90 cm, 120 cm, 240 cm and 500 cm. 



Table 1: Average jam size at different exits of Fig. [6]for the four strategies. 



Door width [m] 


Jam area [m'^] 


LSP 


GSP 


GSQ 


LSQ 


0.90 


53.7 


53.7 


48.8 


49.3 


1.20 


52.8 


52.8 


48.2 


42.5 


2.40 


27.4 


27.4 


31.3 


28.1 


5.00 


13.2 


13.2 


20.3 


18.1 



4. Simulation and Analysis 

The strategies discussed in the previous section are tested on the structure described in Fig. [7| with 
different initial distributions. The simulation area is a simplified model of a section of the ESPRIT arena 
in Diisseldorf, Germany which holds up to 50.000 spectators. The room Rl to R4 are the grandstand with 
dimensions 10 x 20 m^. The tunnels are 2.40 m wide and 5 m long. The exits are 1.20 m wide. The rooms 
R4 and R6 are 10 m wide. The criteria used to evaluate the simulation results are presented in the previous 
section. 



In the first simulation case we are interested in the evacuation of a single block of the stadium. Fig. 7a 
shows the initial homogeneous distribution of this test with 250 pedestrians. The dynamic of the evacuation 
after 60 seconds for the four strategies is presented in Fig. [Sj With the shortest path only two exits 
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are effectively used whereas the quickest path approach takes advantage of fours exits leading to a faster 
evacuation. 



□ DD 
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(a) 250 pedestrians distributed in a single block. 
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(b) 1000 pedestrians distributed in four blocks. 



Figure 7: Investigated facility with initial homogeneous distributions of pedestrians. The total simulation area is 2144 



The results of the evacuation time analysis are summarized in Fig. 9a 1000 runs are performed for 
the evaluation. The dynamic brought by the quickest path leads to a faster evacuation, on average 1 
minute faster. As depicted in Fig. [s] other (less congested or free) exits are used, as will be expected in 
a real evacuation scenario. There is not much difference between the results given by the global and the 
local shortest path, this is due to the shape of the investigated facility. The second value analysed is the 



individual distributions of time spent in jam. The results are presented in Fig. [9b] As expected the results 
are correlated to the evacuation time. With the quickest path, pedestrians spend less time in jam, 30 seconds 
on average. This time is almost tripled using the static strategies. Fig. [lO] shows the evolution of the jam 
size at different exits. Exits without congestions have been left out. The static strategies shows five long 
lasting jams. The quickest path alternatives shows more jams but short lived. The later ones are less critical. 
The total jam size distribution is taken from Fig. 9c The average jam size for the static strategies is 11 
and 9 for the quickest path. 

The second test consists of the simulation of the evacuation of the complete facility with 1000 pedestrians. 



The initial distribution is shown in Fig. 7b A snapshot of the simulation after 60 seconds is shown in Fig. 
11 The results of the jamming time and the evacuation time are summarized in Fig. 12a and Fig. 12b 
respectively. There is no much difference between the mean values of the distributions as in the previous 
case. This is due to the highly congested situation. There is no such "quickest path" in this case and most 
of the pedestrians just have to follow the flow. Those results are confirmed by the time in jam distribution. 
The total jam size distribution is taken from Fig. 12c The mean jam sizes are 46 for the static and 40 



' for the dynamic strategies. 

In the third simulation scenario, the response of the system to a disturbance is simulated. The exits E2 



and E8 (see Fig. 7a) are broken and cannot longer be used (see Fig. 7b). The results of evacuation time 



and time in jam are presented in Fig. |13a| and Fig. 13b respectively. There is now a certain asymmetry 
in the escape route scheme and this is reproduced by the modelling approach. The quickest path lead to 
a faster evacuation and to less time in jam. This behaviour is also the expected one. The total jam size 



distribution is taken from Fig. 13c All mean values are lower than the previous case without disturbance. 
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(a) Local shortest path 
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(b) Local shortest with quickest path 
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(c) Global shortest path (d) Global shortest with quickest path 

Figure 8: Dynamics of the system after 60 seconds for the initial distribution in Fig. [7a| Congestions areas are red. 




(a) Evacuation time distribution (b) Jam time distribution (c) Jam size distribution 

Figure 9: Evacuation time and time in jam distribution for 250 pedestrians. The initial positions are presented in Fig. |7a| 
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(c) Global shortest path 
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(d) Global shortest with quickest path 



Figure 10: Jam size evolution for 250 pedestrians at different exits for the initial dis tribu tion in F ig. \ 7a\ The colours correspond 
to the different exits. Exits without congestions have been left out of the plots. In |(b)| and |(d)| there are more but short-lived 
jams than in|(a)|and|(c)| 
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(a) Local shortest path 
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(c) Global shortest Path 
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(d) Global shortest with quickest path 



Figure 11: Dynamics of the system after 60 seconds for the initial distribution in Fig. |7b| Congestions areas are red. 
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Figure 13: Evacuation time and time in jam distribution for 1000 pedestria ns a fter a system disturbance. Three exits have 
been closed (broken escape route). The initial positions are presented in Fig. |7b| 



5. Summary 

The approach presented in this paper offers the possibihty to assess a given building structure considering 
the manifold of possible route choices. Four strategies have been presented in a combination of quickest 
and shortest path. The quickest path approach, which is based on an observation principle, is not sensitive 
to initial distribution of pedestrians or special topologies like symmetric exits, making it quite general. 
Furthermore it leads to a more realistic dynamics in the evacuation simulation. In addition criteria to 
assess the criticality of an evacuation simulation have been elaborated. We investigated the evacuation 
time distribution, the time in jam distribution and average jam size of pedestrians using a quickest and 
shortest path routing approaches in a graph-based way finding algorithm. The approaches have been tested 
on different scenarios with different complexities involving symmetric, asymmetric or even broken escape 
routes. Similarities between the distributions, their meanings and impacts on different evacuation scenarios 
have been quantitatively analysed. 
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